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Abstract 

Recently there has been an accumulation of experimental evidence in the high temperature superconductors sug- 
gesting the relevance of electron-phonon coupling in these materials. These findings challenge some well-held beliefs 
of what electron-phonon interactions can and cannot do. In this article we review evidence primarily from angle- 
resolved photoemission (ARPES) measurements which point out the importance of electronic coupling to certain 
phonon modes in the cuprates. 
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The physics underlying the peculiar behavior of the 
high temperature superconductors in the normal state 
and the high transition temperatures themselves is far 
from understood. While the importance of magnetic 
interactions has been widely pointed out, a number of 
experimental findings strongly suggest that electron- 
phonon coupling must play a role in these materials. 
These findings include strong phonon renormalizations 
with temperature and/or doping seen in neutron and 
Raman scattering[l,2], the doping dependence of the 
isotope effect of the transition temperature[3], the iso- 
tope dependence of the superfluid density [4], and lastly 
the pressure, layer and material dependence of T c it- 
self. 

In this paper we focus attention on some recent de- 
velopments coming from angle-resolved photoemission 
(ARPES) experiments and theoretical developments of 
anisotropic electron-phonon interactions in general. 

Initially, the attention to bosonic renormaliza- 
tion effects in cuprate superconductors was focussed 
on a "kink" in the electronic dispersion near 50-70 
meV for nodal electrons[5,6,7,8,9] or solely below 
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T c [8,10,ll,12]. This has been further shown to be 
generic to many hole-doped cuprates and is seen both 
in the normal and superconducting states[13]. Re- 
cently, the "kink" phenomenon has been reported for 
electronic states throughout the Brillouin zone (BZ) 
again both above and below T c [13]. Although many- 
body effects reminiscent of strong coupling are known 
to exist below T c in an extended k-space range, the 
observation of kinks in the normal state imposes an 
additional constraint. In contrast to the nodal renor- 
malization which shows little change across T c , data 
in the anti-nodal region reveal a dramatic change in 
the effective coupling through T c [13] , as shown in Fig. 
1 for Bi2Sr 2 Cao.92Yo.o8Cu208+i (Bi-2212). In the su- 
perconducting state, classical Engelsberg-Schrieffer 
signatures of electronic coupling to a bosonic mode are 
seen in the image plot. The image plot shows strong 
" kinks" or breaks in the energy dispersion of the band 
near the (7T, 0) regions of the BZ, and weaker kinks 
near the nodal directions. The shift in the energy 
at which the 40meV mode couples to the electrons 
(4> < 30 deg) compounded by the lack of shift in the 
~ 60 — 70 meV nodal kink (0 > 35 deg) gives rise to a 
fairly uniform ~ 65 — 70 meV kink energy throughout 
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Fig. 1. The image plots in al-a6), and the corresponding MDC-dcrivcd dispersions, arc cuts taken parallel to (0,7r)- (tt, 7r) at the 
locations indicated in the sketched zone at 15K. The image plots in bl-b6), and the corresponding MDC-derivcd dispersions, are 
cuts taken parallel to (0,0) — (7T,7r), at k-spacc locations indicated in the zone at 15K. Kink energies as a function of angle (f> are 
summarized in c) for optimally (94K) doped samples in the normal and superconducting states, d) and e) are spectra taken parallel 
to (0, 7r) — (n, 7r) at the k-space locations indicated for under-doped (85K) and over-doped samples (65K) of Bi-2212, respectively. 



the BZ in the superconducting state. The coupling 
strength, on the other hand, has a strong momentum 
dependence below T c . The relative sharpness of the 
Momentum Distribution Curves (MDC) kink (Fig. 
l(a,b))indicates that the effective coupling strength 
increases from the node to the antinode of the d-wave 
gap. One can also see the increase in coupling in the 
image plots of the raw data towards (-7r, 0) as there 
is a stronger depression in intensity at the mode en- 
ergy. When the coupling gets very strong as in Fig. 
I(al,a2), little dispersion can be tracked, and this is 
the primary indicator of coupling. The kinks occur for 
different dopings, and are clearly seen in underdoped 
as well as overdoped samples, as shown in Figure ld,e, 
respectively. In the deeply overdoped sample (T c 
~ 65K or 5 ~ 22%), Fig. le, the kink energy moves 
to 40-45meV since Ao (~ 10-15meV) becomes much 
smaller. In the underdoped sample (T c ~ 85K), Fig. 
Id, the kink energy remains around 70meV since it has 
a similar gap as the optimally doped sample (Ao ~ 
35-40meV). The signatures of coupling remain strong 
throughout, although they do noticeably increase from 
the overdoped to underdoped sample when comparing 
data taken at the same <j>. 

In the normal state little of these effects can be seen 
without further analysis. As a consequence, kinks in 



the normal state have been missed in prior works. In 
Fig. 2, we have extracted dispersions for three k-space 
cuts in a momentum space region between the nodal 
direction and the Van Hove Singularity (VHS) at (n,0). 
Since we are only concerned with the energy scale, we 
fit the energy dispersion curves (EDCs) phenomeno- 
logically. The usual method to extract dispersions — 
fitting MDC with Lorentzians — is not appropriate 
since the assumed linear approximation of the bare 
band fails towards (n,0) where the band bottom is 
close to Ef[14]. EDC-derived dispersions of three in- 
dependent data sets shown in Fig. 2 consistently reveal 
a ~40meV energy scale in the normal state (al,bl,c), 
evolving into a 70 meV feature in the superconducting 
state (a2,b2), as summarized in Figure 2d. In the su- 
perconducting state where the coupling is strongest, 
the peak position in both EDC and MDC dispersions 
asymptotically approach the characteristic energy 
defined by the bosonic mode (Fig. 2a2, Fig. 2b2). 

The kink energy shift is due to the opening of a su- 
perconducting gap which shifts the energy at which the 
electronic states couple to the bosonic mode. Here, we 
observe a kink shift of ~25-30meV, close to the maxi- 
mum gap energy, Ao. We summarize the temperature 
dependence of the energy at which we see a bosonic 
mode couple to the electronic states in the anti-nodal 
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Fig. 2. EDC(al, bl, c) derived dispersions in the normal state 
(107K and 115K). <fi and the cut-direction are noted in the 
insets. The red dots are the data; the fit to the curve (dashes, 
black) below the 40meV line is a guide to the eye. a2) and b2) 
are MDC derived dispersions at the same location and direction 
as in al) and bl), but in the superconducting state (15K). In 
a2) and b2) we also plot the peak(I) and hump positions(II) of 
the EDCs for comparison. The inset of a2) shows the expected 
behavior of a Bogoliubov type gap opening. The s-likc shape 
below the gap energy is an artifact of how the MDC handles the 
back-bend of the Bogoliubov quasiparticlc. d) kink positions 
as a function of (p in the anti-nodal region. 



region in Figs. 1(c) and 2(d): the kink energies are at 
~40meV above T c near the anti-nodal region and in- 
crease to ~70meV below T c . Because the band mini- 
mum is too close to Ef , the normal state ~40meV kink 
cannot easily be seen below <j> ~ 20°. 

By comparing the different cut directions, one can 
clearly see that the coupling is extended in the Bril- 
louin zone and has a similar energy scale throughout, 
near ~ 70meV. Moreover, the signatures of coupling 
increase significantly toward (it, 0) or smaller 0. The 
clear minimum in spectral weight seen near 70meV in 
Fig. 2al) and 2a2) indicates strong mixing of the elec- 
tronic states with the bosonic mode where the two bare 
dispersions coincide in energy. 

The above described experimental observations are 
difficult to reconcile in the spin resonance mode sce- 
nario. The kink is clearly observed at 40meV in the 
normal state at optimal doping while the 41meV spin 
resonance mode exists only below T c [15]. The kink is 
sharp in the superconducting state of a deeply over- 
doped Bi2212 sample (Fig. le, consistent with data of 
[12]) while no spin resonance mode has been reported, 
or is expected to exist at this doping. In addition the 
kink seen in under-doped Bi2212 (Fig. Id) is just as 
sharp as in the optimally doped case, while the neutron 
resonance peak is much broader[15]. Lastly, it is diffi- 



cult to account for the strong kink effect seen through- 
out the BZ from the spin mode since the mode's spec- 
tral weight is only 2% [16] as it may only cause a strong 
enough kink effect if the interaction with electrons is 
highly concentrated in k-space [17] . 

On the other hand, it has been argued that the 
bosonic renormalization effect may be reinterpreted as 
a result of a coupling to the half-breathing in-plane Cu- 
O bond-stretching phonon for nodal directions[6,7] and 
the out-of-plane out-of-phase O buckling Bi g phonon 
for anti-nodal directions[13]. This re-interpretation has 
the clear advantage over the spin resonance in that it 
naturally explains why the band renormalization effect 
is observed under three circumstances: 1) in materials 
where no spin mode has been detected, 2) in the nor- 
mal state, and 3)in the deeply overdoped region where 
the spin mode is neither expected nor observed. How- 
ever this interpretation also requires that the electron- 
phonon interaction be highly anisotropic and its im- 
pact on the electrons be strongly enhanced in the su- 
perconducting state. This is something one does not 
expect a priori. 

We now show that the anisotropy of the electron- 
phonon interaction of these phonons, when formulated 
within the same framework as that used to explain 
phonon lineshape changes observed via Raman and 
neutron measurements[l,2,18,19], explain the ob- 
served momentum and temperature dependence of 
the data, leading to unified understanding of band 
renormalizations. The electron-phonon interactions 
are very anisotropic as a consequence of the following 
four properties of the cuprates: 1) the symmetry of 
the phonon polarizations and electronic orbitals in- 
volved leading to highly anisotropic electron-phonon 
coupling matrix elements g(k,q); 2) the kinematic 
constraint related to the anisotropy of the electronic 
band structure and the van Hove singularity (VHS); 
3) the d-wave superconducting gap, and 4) the near 
degeneracy of energy scales of the B\ g phonon, the 
superconducting gap, and the VHS. 

We consider a tight-binding three-band Hamiltonian 
modified by Bi g buckling vibrations as considered in 
Ref. [18] and in-plane breathing vibrations as consid- 
ered in Ref. [20]: 

n,a 

n,<5,<j 

+ X Ps[t - Q«5d P w?(n)](&J 1 ^a n ,«,<7 + h.c.) 
n,s,a 

+t' ^ Ps,s>i a ii,s a n,8',<T + h.c) 
n,s,s' 



(1) 



with Cu-0 hopping amplitude t, O-O amplitude t' , and 
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Fig. 3. Plots of the elcctron-phonon coupling | p(k, q) | for 
initial k and scattered — k — q states on the Fermi sur- 
face for the buckling mode (left panels) and breathing mode 
(right panels) for initial fcrmion k at an anti-nodal (top pan- 
els) and nodal (bottom panels) point on the Fermi surface, 
as indicated by the arrows. The red/blue color indicates the 
maximum/minimum of the cl-ph coupling vertex in the BZ for 
each phonon. 

e^,p denoting the Cu and O site energies, respectively. 
Here a a - XlV , a^ a=x y (b,b^) annihilates, creates an elec- 
tron in the a (Cu) orbital, respectively. The overlap 
factors Qs,Ps, and P' s S i are given in our notation as 
Q±x = Q± y = ±l;P±a = ~P±y ~ ilj P±x,±y — 1 — 



~f%x,±y F° r the half-breathing mode we consider only 
couplings gd p = qot (with 1/qo a length scale) arising 
from modulation of the covalent hopping amplitude 
£[20,21]. For the Bi g mode a coupling to linear order 
in the atomic displacements arises from a local c-axis 
oriented crystal field E z which breaks the mirror plane 
symmetry of the Cu-0 plane[18]. Such a field may be 
due to static buckling or different valences of ions on 
either side of the Cu-0 plane, orthorhombic tilts of the 
octahedra, or may be generated dynamically. 

The amplitudes 4> of the tight-binding wavefunctions 
forming the anti-bonding band with energy dispersion 
e(k) arising from Eq. 1 include the symmetry of the p, d 
wavefunctions of the Cu-O plane. The specific trans- 
formation is constructed from the projected part of 
the tight-binding wavefunctions 6^ CT = 0(,(k)di <; n and 
a a ,k,cr = a (k)d k ct , 



4>x,y (k) 



1 



JV(k) 



JV(k) 



e 2 (k)-t' 2 (k)], (2) 



[e(k)t x , y (k)-t'(k)t y , x (k)], (3) 



with iV 2 (k) = [e 2 (k) - t' 2 (k)} 2 + [e(k)t«(k) - 
t'(k)t y (k)] 2 + [e(k)t a (k) - t'(k)t x (k)] 2 , t a (k) = 
2i sin(fc a a/2), and i'(k) = W sm(k x a/2) sm(k y a/2). 
In particular if k resides on the Fermi surface, then 



>f S (k) 



f(k) 



^t>(ky + t x (ky + t y (ky 



(4) 



±isgn[t'Ck)}ty !X (X) 



(5) 





^1 Buckling 

^ BO = 0.23 

H Brealhing 







1° 

0.52 
1.0 
1.6 

' 2.6 

l° 

0.060 

0.12 
0.18 
( 0.24 
' 0.30 



Fig. 4. Plots of the electron-phonon coupling Ak in the first 
quadrant of the BZ for the buckling mode (right top panel) 
and breathing mode (right bottom panel). The color scale is 
shown on the right for each phonon. The left panel shows 
energy contours for the band structure uscd[23]. 

The full anisotropy of the electron-phonon cou- 
pling comes from the projected wavefunctions in 
combination with the phonon eigenvectors: H e i- P h = 
^£k,q, CT , M 3M(k,q)c^ CT C k+a >^ q + a M ,_q]. Ne- 
glecting the motion of the heavier Cu atoms compared 
to O in the phonon eigenvectors the form of the re- 
spective couplings can be compactly written as 



9B lg (k, q) = eE z 



h 



(6) 



4M M(q)fi fll8 
x {<f> x (k)<f) x (\<i) cos(q y a/2) — (f>y(k)cf)y(]s.') cos(q x a/2)} 

(7) 



gv(k,q) = g dp 



h 



2Mo£lbr 



E 



x {<j) b (k')<t> a (k) cos(k' a a/2) — 0(,(k)0 a (k') cos(k a a/2)} 

with M(q) = [cos 2 (ga;a/2) + cos 2 (gj,a/2)]/2 and 
k' = k-q. Here we have neglected an overall irrele- 
vant phase factor. The anisotropy can be shown at 
small momentum transfers q where (?s lg (k, q = 0) ~ 
cos(k x a) — cos(k y a), while gtr ~ sin(qa) for any k. 
Also for large momentum transfer q = (it /a, it /a), 
9B lg vanishes for all k while gi, r has its maximum for 
k located on the nodal points of the Fermi surface. 
While this g-dependence has been the focus before in 
context with a d x 2_ y 2 pairing mechanism [22], the de- 
pendence on fermionic wavevector k has usually been 
overlooked. This strong momentum dependence of the 
Big coupling implies that resistivity measurements 
would be sensitive only to the weaker part of the 
electron-phonon coupling, in agreement with experi- 
ments. It is exactly this fermionic dependence which 
we feel is crucially needed to interpret ARPES data. 

Fig 3 plots | g(k, q) | 2 as a function of phonon 
scattered momentum q connecting initial k and final 
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Fig. 5. Image plots of the calculated spectral functions in the normal (al,bl,cl) and superconducting (a2,b2,c2) states compared 
to the spectral functions in the normal (a3,b3,c3) and superconducting (a4,b4,c4) states measured in Bi2Sr2Cao.92Yo.o8Cu20s+5 
(Bi-2212)[13] for momentum cuts a,b,c shown in the right-most panel and in Figure 2. The same color scale is used for the 
normal/superconducting pairs within each cut, but the scaling for the data and the calculation are separate. The red markers 
indicate 70meV in the superconducting state. 



k' = k — q fermion states both on the Fermi surface. 
Here we have adopted a five-parameter fit to the band 
structure e(k) used previously for optimally doped Bi- 
2212[23]. For an anti-nodal initial fermion kAJV, the 
coupling to the Bi g phonon is largest for q — and 2kj? 
scattering to an anti-nodal final state. The coupling to 
the breathing mode is much weaker overall, and is espe- 
cially weak for small q and mostly connects to fermion 
final states at larger q. For a nodal fermion initial state 
kjv, the Bi g coupling is suppressed for small q and has 
a weak maximum towards anti-nodal final states, while 
the breathing mode attains its largest value for scat- 
tering to other nodal states, particularly to state — kjv- 
In Fig. 4 we plot Ak for each k point in the BZ in the 
normal state, where we have defined [24]: 

* = 5^ £[/(<«) -/(*<+«i)>] 

q 

x I 5 (k, q) | 2 <5(e(k + q) - e(k) - Klo), (8) 

with Nf = jf X^k ^( e ( k )) the density of states per spin 
at the Fermi level. As a result of the anisotropy of the 
coupling constant, large enhancements of the coupling 
near the anti-node for the B\ g phonon, and nodal di- 
rections for the breathing phonon mode are observed, 
respectively, which are an order of magnitude larger 
than A = i ^ k Ak- Here we have used the value of the 
electric field eE z = 1.85eV/A consistent with Raman 
scattering measurements of the B\ g phonon and buck- 
ling of the Cu-O plane in YBaCuO, and have taken 
qot — 2^/2eE z as a representative value for the cou- 
pling to the breathing mode. These values lead to net 
couplings similar to that obtained from LDA calcula- 
tions for an infinite layer compound CaCu02 which 
develops a static dimple [24]. 

The Nambu-Eliashberg electron-phonon self- 
energy determines bosonic features in the spectral 



function[25]: E(k,io; ra ) = iu) m (l — Z(k, iw m ))fo + 
X(k,iw m )f3 + $(k,ia; n ,)fi with 

S(k,iw m ) = — ffM( k >P- k )3M(P> k -P) 
xD°(k-p, iw m - iuj n )T 3 G°(p, iuj m )r 3 , (9) 

with G°( V ,iu m ) = '" m ;L + ^)?-^(p ( ) p)fl and s2 ( k ) = 

e 2 (k) + A 2 (k), with Ti =u ..-3 are Pauli matrices. We 
take A(k) = Ao[cos(fc x a) — cos(k y a)]/2 with Ao = 
35meV. The bare phonon propagators are taken as 
D°(q,iQ v ) = | [ mJl^ - • Here we take 

dispersionless optical phonons h£lB lg ,br = 36, 70meV, 
respectively. 

The spectral function A(k,uj) = — ~ImGu(k, iu) — > 
to + iS) is determined from Eq. 9 with G(k,iu) m ) = 
[foiu; m Z(k, iujm) — (e(k)+x(k, itj m ))f 3 — <l>(k, io; m )fi] _:1 
and the coupling constants given by Eq. 6-7. Here we 
consider three different cuts as indicated in Fig. 4: 
along the nodal direction (k x = k y , denoted by c), 
and two cuts parallel to the BZ face (k y = 0.75n/a, 
b) and (k y = O.Gin/a, a) where bilayer effects are not 
severe [13]. 

Our results for the three cut directions are shown 
in Fig 5 for the normal (T = 110K) and supercon- 
ducting (T = 10K) states. Very weak kink features 
are seen in the normal state, with a smeared onset of 
broadening at the phonon energy due to the measure- 
ment temperature. A much more pronounced kink at 
70 meV in the superconducting state is observed for 
all the cuts along with a sharp onset of broadening, in 
excellent agreement with the data[13], shown also in 
Fig. 5. In the superconducting state, the kink sharp- 
ens due to the singularity in the superconducting den- 
sity of states from the anti-nodal regions. Thus cut a 
exhibits the strongest signatures of coupling, showing 
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an Einstein like break up into a band that follows the 
phonon dispersion, and one that follows the electronic 
one. In cut b, a weaker signature of coupling manifests 
itself in an s— like renormalization of the bare band that 
traces the real part of the phonon self-energy, while in 
cut c, where the coupling to the Bi g phonon is weak- 
est, the electron-phonon coupling manifests itself as a 
change in the velocity of the band. The agreement is 
very good, even for the large change in the electron- 
phonon coupling seen between cuts a and b, separated 
by only l/10th of the BZ. 

We now discuss the reason for the strong anisotropy 
as a result of a concurrence of symmetries, band struc- 
ture, d— wave energy gap, and energy scales in the 
electron-phonon problem. In the normal state, weak 
kinks are observed at both phonon frequencies, but 
more pronounced at 70 meV for cut c and at 36 meV 
for cut a, as expected from the anisotropic couplings 
shown in Fig. 3. However this is also a consequence of 
the energy scales of the phonon modes and the bottom 
of electronic band along the cut direction as shown in 
Fig. 4. Further towards the anti-node, another inter- 
esting effect occurs. As one considers cuts closer to the 
BZ axis, the bottom of the band along the cut rises in 
energy and the breathing phonon mode lies below the 
bottom of the band for (k x = 0, k y > 0.827r/a) and so 
the kink feature for this mode disappears. This is the 
usual effect of a Fano redistribution of spectral weight 
when a discrete excitation lies either within or outside 
of the band continuum [26]. This thus can be used to 
open "windows" to examine the coupling of electrons 
to discrete modes in general. Taken together, this nat- 
urally explains why the normal state nodal kink is near 
70 meV[5,6,7,8,9] while the anti-nodal kink is near 40 
meV[10,ll, 12,13]. 

In the superconducting state the bare band renor- 
malizes downward in energy due to the opening of the 
gap and gives the strongest effect for the anti-nodal re- 
gion where the gap is of the order of the saddle point en- 
ergy. The Bi g coupling becomes dramatically enhanced 
by the opening of the gap at a frequency close to the 
frequency of the phonon itself, and leads to the dra- 
matic renormalization of the phonon observed in Ra- 
man and neutron measurements [1]. Yet the breathing 
coupling does not vary dramatically as the large gap 
region is not weighted as heavily by the coupling con- 
stant. As a consequence, the Bi g mode dominates and 
clear kinks in the data are revealed at around 70 meV 
- the -Bi 9 energy plus the gap. Indeed a much weaker 
kink at 105 meV from the breathing phonon shows at 
higher energies for the nodal cut but this kink becomes 
weaker away from nodal directions as the coupling and 
the band bottom along the cut moves below the phonon 
frequency so that the "energy window" closes. 

Finally we plot in Fig. 6 the superconducting density 
of states (DOS) ~l/nImGn(k,ui) calculated with the 
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Fig. 6. Superconducting density of states calculated with the 
electron-phonon self energy for three different values of the 
electric field eE z — (black), 1.85 (blue), and 3.2 (red), in units 
of eV/A, respectively. The yellow dashed lines indicate the 
energies ^isig + Ao and + Ao- 

electron-phonon contributions from both the breath- 
ing and Bi g phonons to the self energy determined via 
Eq. 9 for three different values of the Bi g coupling con- 
stant parametrized by the crystal field eE z . We have 
fixed the relative breathing coupling at qot = 2^/2eE z 
as for ARPES. For these values of the coupling con- 
stants the structure in the DOS is determined largely 
by the anti-nodal couplings to the Bi g phonon. For 
the value of the coupling which fits the ARPES data 
well, the peak-dip-hump structure of the DOS (at val- 
ues E = Ao, f2si 9 + Ao, and Q,Bi g + E v h, respectively) 
compares very well with the DOS determined via scan- 
ning tunnelling microscopy in Bi-2212[27]. 

In summary, contrary to usual opinion of the role 
of anisotropy in electron-phonon coupling, the inter- 
play of specific coupling mechanisms of electrons to the 
buckling and breathing phonons give a natural inter- 
pretation to the bosonic renormalization effects seen 
in ARPES in both the normal and superconducting 
states, and provides a framework to understand renor- 
malizations as a function of doping. 
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